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Abstract 

We study the equilibration of the right-helicity states of light Dirac neutrinos 
in the early universe by solving the momentum dependent Boltzmann equa- 
tions numerically. We show that the main effect is due to electroweak gauge 
boson poles, which enhance thermalization rates by some three orders of mag- 
nitude. The right-helicity states of tau neutrinos will be brought in equilibrium 
independently of their initial distribution at a temperature above the poles if 
rriu, > 10 keV. 
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1 Introduction 



Primordial nucleosynthesis is a remarkable probe of neutrino properties . Although 
recently the increasing accuracy of the cosmological data, such as observations related 
to the primordial abundances of helium, deuterium and the other light elements, has 
emphasized systematic errors inherent to primordial nucleosynthesis analysis, there 
remains a great potential for constraining neutrino physics via cosmological observa- 
tions. To some extent primordial nucleosynthesis could be sensitive even to the Dirac 
vs. Majorana nature of neutrinos, because in the Dirac case the small relic abundance 
of the inert right-handed component of Dirac neutrino would also contribute at nucle- 
osynthesis. Of course, presently one cannot hope to differentiate between the Dirac 
and Majorana nature of neutrinos on cosmological grounds, but in principle this is 
an interesting problem. In practise, because of the smallness of the electron neutrino 
mass, from this point of view only the right-handed components of and Uj. can have 
interesting relic abundances. 

Naively, the cosmological density of light right-handed neutrinos (or rather the 
right-helicity states of light neutrinos, z/_|_) is expected to be very small. If i/r's ever were 
in equilibrium, they decoupled very early because of their low capability of interacting. 
Assuming this took place well above the electroweak phase transition temperature 
Tew, at the onset of primordial nucleosynthesis the contribution of the right-handed 
component of a Dirac neutrino is given by pR ~ [g*{l MeV)/(7*(TEw)]^''^PL — 0.044pi, 
where g*{T) is the effective number of degrees of freedom in thermal equilibrium in 
the temperature T and pi is the equilibrium energy density of left-handed neutrino. 
This assumes that at high temperature we may consider the gas of quarks, leptons, 
and gauge bosons nearly ideal, which may not be true. A recent lattice simulation 
of the QCD energy density above the critical temperature has revealed that the actual 
energy density is some 15% smaller than expected, which might signify the existence 
of a condensate at high T. One also assumes that there is no significant entropy 
production either at the electroweak or QCD phase transitions. Lattice simulations 
seem to indicate that this is true for QCD, and the latent heat in the electroweak 
phase transition is also known to be small . 

The right-helicity states of Dirac neutrinos are not completely inert in the Standard 
Model but can be produced (and destroyed) in spin-flip transitions induced by 
the Dirac mass § or the neutrino magnetic moment. If the neutrino mass is 
large enough, i/+ would be produced in collisions below the QCD phase transition 
temperature, resulting effectively in an additional neutrino species at nucleosynthesis. 
An important source of z/+'s are also the non-equilibrium neutrino scatterings and 
decays of pions, as was pointed out in 0. This gives rise to the bound m,,^ < 130 keV 
and rriiy^ < 150 keV pj, using Tqcd = 100 MeV and assuming that nucleosynthesis 



allows less than 0.3 extra neutrino families. This is very restrictive bound for tau 
neutrinos since, basing on primordial nucleosynthesis arguments, there seems to be 
no window of opportunity for a sufficiently stable (tj, > 10^ sec) tau neutrino in the 
MeV region because of the production of non-equilibrium electron neutrinos in i/^-z/^ 
annihilations 0. 

The actual cosmological density of the right-helicity neutrinos depends not only 
on the production rate near the QCD phase transition, but also on whether the right- 
helicity neutrinos had a chance to equilibrate at some point during the course of the 
evolution of the universe. In this paper we wish to point out that at T < 100 GeV 
there is an enhancement of the z/+ production rate due to the electroweak gauge boson 
poles.Q As a consequence there is a temperature range in which the production rate 
can exceed the expansion rate of the universe and the right-helicity neutrinos may be 
brought into thermal equilibrium with other light particles. This will take place if 
the neutrino mass is sufficiently large. We show that for the tau neutrino the mass 
limit is about 10 keV, and of the same order of magnitude for the muon neutrino. At 
the temperature T = 1 MeV the energy density of right-helicity tau neutrinos with 
a mass of 10 keV is found to be about 6% of the energy density of an ordinary left- 
helicity neutrinos, i.e. their contribution to the effective number of neutrino species is 
AA^,^ ~ 0.06. The contribution of right-helicity neutrinos with a smaller mass depends, 
apart from the mass, on their initial energy density above the electroweak scale. 

The paper is organized as follows. In Section 2 we will list all tree-level spin-flip 
reactions where right-helicity neutrinos are produced and discuss their relative im- 
portance. We will also demonstrate the ly-pole effect by considering the reactions 
ud —>■ z/+r+ and t^u —>■ v^d as an example for the right-helicity tau neutrino pro- 
duction. In Section 3 we consider the evolution of right-handed neutrino density and 
describe the numerical method used in solving the Boltzman equation. The results 
and a discussion is presented in Section 4. 

2 Right-helicity neutrino production 
2.1 Processes 

We shall consider the production of right-helicity neutrinos in the Standard Model 
without assuming any new interactions except the Yukawa interactions of the right- 
handed neutrino with the scalar doublet. The Yukawa coupling is the origin of Dirac 
neutrino mass and provides a spin-flip operator responsible for the interactions of 
the right-helicity neutrinos, the probability of the spin flip being proportional to the 

^ The question of the cosmic abundance of the right-handed component of Dirac neutrinos was 
previously studied in There the pole effect was not considered. 



neutrino mass squared. Another source for the spin-flip is neutrino magnetic moment 
[0], which appears in the Standard Model first at one- loop level and is, therefore, too 
small to be of any significance for our considerations. 

We shall consider light neutrinos with mass less than 0.2 MeV at temperatures 
between 1 MeV and 100 GeV. Since each right-helicity neutrino introduces in the 
matrix element a factor my/\p^\, we may safely neglect processes in which more than 
one right-helicity neutrino is involved. All relevant processes are listed (up to crossing) 
in Table 1. 

There are 68 purely fermionic 2 — 2 processes in which a right-helicity muon or 
tau neutrino can be produced. In addition, a right-helicity tau neutrino can also be 
produced in 11 lepton and quark three-body decays, and the muon neutrino in another 
set of 11 three-body decays. Since we are especially interested in interactions occuring 
around the poles of weak gauge bosons, we have to consider also processes involving 
W^, Z and H. There are 16 such processes. Finally, there are 3 two-body decays 
of W^, Z and H bosons which are capable producing right-helicity muon and tau 
neutrinos. However, as we will show below, the processes involving gauge or Higgs 
bosons can be neglected in comparision with the purely fermionic processes. 

In what follows, for definiteness, we will consider only processes including right- 
helicity tau neutrinos. Since we are mainly interested on processes at temperatures 
above the muon and tau lepton masses, we expect that the results below are roughly 
valid also for muon neutrinos. 



2.2 Production rates 

Let us consider a 2 ^ 2 scattering a + h v^ + d where one of the final state particles 
is a right-helicity neutrino. To estimate the relative importance of various processes 
we approximate the thermally averaged production rate per one z/+ by 

r+ = I dIiadIibdIi+dnd{2Tx f6'^''\pa+Pb-p+-Pd)S\Mab^+d? 

Ti_|_ J 

^fmi-fDii-fD, (1) 

where n™ is the the equilibrium number density of the right-handed neutrinos, dUi = 
d^Pi/{{27r)^2Ei), S is the symmetry factor taking into account identical particles in 
the initial and/or final states, and /™ are Fermi-Dirac distribution functions. 

At high T the rate Eq. (|I]) is infrared sensitive to the thermal corrections in the 
propagators. In general, the structure of the gauge boson propagators at finite T 
is complicated because Lorentz symmetry is lost, which results in separate transverse 
and longitudinal self energies IIt{^, k) and k). In practise, however, the leading 

thermal effect arises from small momenta, so that in most cases it is an excellent 
approximation just to modify the propagators by introducing a Debye mass M^(T) = 



UL{uj,k = 0), which we approximate by Mf(T) ~ + 0.1 (i = W,Z). (This 
modification is necessary in t-channel propagators only, since s-channel propagators 
are not infrared sensitive.) In this approximation external particles or interaction 
vertices do not receive thermal corrections. Thus the dispersion relations in the thermal 
distributions in Eq. (0) remain unchanged. Admittedly, this is a simplistic approach, 
but for our purposes, and for the desired accuracies, this should be sufficient. In fact, 
in the region of interest the effects due to thermal masses turn out to be very small. 

In addition to thermal corrections, we must account for the imaginary parts of 
the gauge boson propagators, or the widths. This is particularly important for the 
s-channel. Thermal corrections will generate additional imaginary parts both in the 
s-channel and t-channel, but at T < 100 GeV they may safely be neglected. 

A technical detail worth pointing out is that for a fixed helicity, the spin-fiip matrix 
elements are not Lorentz-invariant since the direction of the spin picks out a preferred 
frame of reference, as was emphasized in Indeed, a Lorentz boost changes the 
helicity of the particle, so that sometimes a fixed-helicity reaction forbidden in the 
CM may actually take place in another frame. This means that it is not sufficient to 
compute matrix elements just in e.g. the CM-frame, but instead one should a use a 
general frame. 

The main purpose of this paper is to show that interactions at the weak boson pole 
may bring the right-helicity neutrinos to thermal equilibrium. To demonstrate this 
effect, let us consider the the t-channel reaction t~u — > u+d and its crossed s-channel 
reaction ud , where z/+ is a right helicity tau neutrino. The matrix element 

for the t-channel process reads 

M = ^VudRw{q'){-9,u + ^) u,Y{l - l,)uru,r{l - l6)uu , (2) 
V 2 M^/ 

where q = Pr — Pu, Vud is the appropriate CKM matrix element, and as in s-channel, 
the propagator gives rise to a term of the form 

q^ - M^{T) + iTwMw 

As mentioned above, the width T\y is important only for the s-channel. Note the 
T-dependence in the numerator. The matrix element squared summed over spins of 
T~ and quarks is then 

\M\'' = Gl\Vud?\Rw{q^)\\Tgg + T,, + T,,) , (4) 



where 



Tgg = 64:{Pr ■ Pu){Ku ■ Pd) , 

Tgq = ^ { [mliPu ■ Pd) - mliPv ■ Pu)] ■ Pr) 
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X [[ml + ml] {pu ■ Pd) - 2mlml] . 



(5) 



Here K = p — ms , with the spin four- vector s for particles with a definite hehcity 
h given by 
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In the ultra-relativistic hmit K''^ can be approximated as 



2p^ for h = -1 
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Accordingly the quantities Tgg, Tgg and Tgg in Eq. have in the case of z/+ production 
typical sizes given by 
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For the production of one has 
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Tgq ~ 



K/M, 



and Tgq ~ 



{niqmr / M'y/Y'T'^ . For both z/+ and z/_ at T < My/ and in the case of quarks other 
than the top, the terms Tgq and Tqq can be neglected in comparison with Tgg. The 
rates for z/+ are suppressed by a factor of the order of ml/T"^ compared with those 
for z/_. However, at higher temperatures, T ^ Myy, the Tgq term starts to dominate 
the production rate of because of terms not proportional to K,^. Consequently the 
production rates of z/+ are in this case suppressed by a constant factor of the order of 
m^ml/M^. 

The s-channel matrix element can easily be obtained from Eq. (^) by crossing. 
The thermal rates for both s- and t-channel processes can then be found from Eq. (|l|). 
A numerical integration results in the curves displayed in Fig. 1. Here we have for 
simplicity ignored the small finite temperature effects. It can be see from the figure 
that the effect of the pole is spread over a relatively large temperature range. This is 
a consequence of thermal averaging. Nevertheless, the enhancement in the s-channel 
is apparent. The results presented in Fig. 1 are for the tau neutrino of mass 20 keV, 
but they are equally valid for the muon neutrino. In this case the thermally averaged 



production rate r+ of the s-channel process is seen to exceed in a certain temperature 
interval the expansion rate of the universe, given by the Hubble parameter 



Hence, with this reaction alone a 20 keV right-helicity tau (and muon) neutrino would 
be brought into thermal equilibrium while universe cools through this stage. The 
complete analysis described in the next section, which is based on solving the Boltz- 
mann equation with all the relevant processes included, confirms this expectation. The 
pole effects are indeed important for an estimate of the relic density of right-helicity 
neutrinos. 

In addition to fermionic 2^2 scattering processes, at high T a potential source 
for z/+'s are 2^2 processes that involve the gauge or the Higgs bosons in the final 
or initial state (we dub such processes "bosonic"). Let us consider as an example 
the process r^7 u^W', where the photon couples either to the charged lepton 
(Compton scattering) or to the W-boson in a three-boson vertex. For comparision, 
the thermally averaged rate of this process is also displayed in Fig. 1. One can see 
that it can be neglected even at temperatures around the pole, because there the 
production rate of the s-channel purely fermionic process is about three orders of 
magnitude higher. 

We have not considered all possible bosonic processes. It is however very plausible 
that generically among the 2 — > 2 scattering processes only the purely fermionic pro- 
cesses are important, and among these, s-channel dominates over t-channel because of 
the pole in the s-channel. In what follows, we will always disregard bosonic processes. 

The importance of decays is less straightforward to discern. Numerical inspection 
reveals that the two-body decays hsted in the Table 1 can be neglected: their contribu- 
tion to the total rate is ~ O{10~^^) at T > 10 GeV, at lower temperatures their contri- 
bution vanishes exponentially. Three-body decays are more important. We find that 
while the total contribution from three-body decays is negligible at higher momenta, 
Ip+I/T" > 3, it is as large as few per cents for \p^\/T ~ 3 at T < 1 GeV, increasing 
up to ~ 30% for very small momenta. For such small momenta also production 
by t — > fer^i/+, which has its maximal contribution (~ 10%) around T ~ 30 GeV, is 
important. 

3 Evolution of the right-helicity neutrino density 
3.1 The Boltzmann equation 

Our goal is to estimate the relic abundance of the right-helicity neutrinos func- 
tion of time. Of particular interest is the contribution of the right-helicity neutrinos 




(10) 



to the effective number of neutrinos, AA^^,, at the onset of primeval nucleosynthesis. 
To find AA^,^ one has to determine first the evolution of the phase space distribution 
of the right-helicity neutrinos t)Q Then using a primordial nucleosynthesis 

code one computes the increase in the primordial ^He abundance Al^ resulting from 
the non-zero energy densities of z/+ and z7_ at the nucleosynthesis time. Since "^He 
abundance is a monotonic function of the total energy density of the universe, one 
can translate Al^ back into AA^",^, the effective change of the energy density in units 
of that of massless two-component neutrinos. Without using the NS code one can 
approximate AA^,^ as 

AAT, ~ p^^^^J—T^^^ , (11) 

where the subscript n ^ p refers to the freeze-out of the reactions which transmute 
protons and neutrons into each other. The justification of this approximation comes 
from the fact that Yp is predominantly determined by the density of neutrons just after 
the n ^ p freeze-out. 

Let us briefly describe the method we have applied for solving the evolution of the 
right-helicity neutrino density. The evolution of the distribution /+(|p_,_|, t) is governed 
by the Boltzmann equation 

with the initial condition 

MlP+lto) = fUM . (13) 

Instead of time t, we will consider the cosmic scale factor R as the independent variable 



describing the evolution. Eqs. (|1^, [T^) then transform into 




\p^\Jn]U = i7{^-^] (14) 




coll 

and 

U{MR{to)) = f4M . (15) 

Since R is defined up to a multiplictive constant (only the ratio of R's at two different 
times has physical meaning), one may choose -R(to) = 1- 

For solving Eqs. ( |T4|, |T5|) it is important that the first-order PDE ([l^) may be 
transformed to an ODE by a suitable transformation of the variables \p_^ \ and R. We 
use ^ 

\P+\ ^ P+ = \P+\-B- ^ (16) 

Kq 



^The arguments of the distribution function /_|_(|p+|,i) indicate that we work in the FRW, i.e., 
spatiaUy isotropic and homogeneous cosmology. 



where Rq is arbitrary but fixed. Then Eq. ( p^ may be written as 

(IT) 

where, according to (p^), /+(p+,-R) := = p+Ro/ R, R) and a similar relation 

holds between the collision terms {df+/dR)coii and {df+/dR)co\\- Eq- (P^ describes 
the evolution of f^{p^, R) with an increasing R for a fixed value of the parameter p^. 
For noninteracting zz+'s, Eq. (|T7p is solved by /+ = const, which may be translated 
into the familiar evolution (caused by momentum redshift) of the distribution function 
of the freely expanding gas of particles, 

U{\p^lR) = flip^ = \p^\R/Ro). (18) 

This equation describes the compression of the distribution function along the mo- 
mentum axis as R increases. For ultrarelativistic particles this compression preserves 
the shape of the distribution function. Interactions will in general introduce distor- 
tions in the distribution function because of the momentum-dependent strength of the 
coupling of particles with the ambient matter. 

As has been already mentioned in Sec. 2.1, we are considering ultrarelativistic 
neutrinos and therefore are justified to neglect all processes involving more than one 
right-helicity neutrino. Then, as will be shown in the next subsection, neglecting for a 
moment a slight /^.-dependence of the Hubble parameter H and the R—T relationship, 
used in calculation of both H and the collision term, the r.h.s. of the Eq. (|l^) is 
linear in the distribution function of right-helicity neutrinos /+, and consequently (p!7|) 
represents a set of uncoupled differential equations for f+{p^,R), one independent 
equation for each value of the parameter p+.^ In our computations we followed the 
evolution of the distribution function at 30 different values (~ bins) of p+ placed at 
equal distances in the interval < p+/(100 GeV) < 10. We took into account the weak 
dependence of H and T{R) on /+ by solving Eq. (^) iteratively, that is, we evaluated 
these quantities using the value of /+ obtained in the previous iteration. Because of 
the smallness of the correction to H and T{R) from /+ the solution converges rapidly 
and we had to use only 3 iterations. 



3.2 The collision term 

Considering only purely fermionic 2^2 and 1 — > 3 processes, we can write the 
collision term on the r.h.s. of ([T7|) in the form 

= {C2^2 + C^^3){1 - /+) - (CU + C[^3)h , (19) 
coll 

^ In the general case when one takes into account also processes with more than one right-helicitity 
neutrinos, the r.h.s. is no more linear in /+ and one arrives at the set of coupled ODEs for /+: the 
evolution of f+{p+, R) and R), for 7^ are not independent. 




where the coefficients Cj represent production and Cj destruction of u^, with 1 = 2 
2, 1 3. The exphcit expressions for the quantities 6*2^2 and Ci^s are 



C' 



Ci^3(|p+|, i?) = E ^ / dUfdU,dn,{2n)'6^'\pj - p, - p+ - pn) 

X S\Mf^,^u?ff''{l - - /D , (20) 



dec 



where, as before, dlli = d^pi/{{27T)^2Ei) and is a symmetry factor taking into account 
identical particles in the initial and/or final states. Note that these expressions do not 
include the factors Qi representing the number of spin degrees of freedom. According 
to our convention these factors are included already into the matrix elements squared 
since we sum over polarization states of all particles except the right-helicity neutrino 
under consideration. We have also used the well justified assumption that charged 
leptons, quarks and left-helicity neutrinos have thermal FD distributions. 
The expressions for Cj can be obtained from (^) by the replacement 

(21) 

Using (^OD and (^) one easily finds that the primed and unprimed coefficients are 
related to each other by 

C'j = e^+/^Ci , (22) 

where we have assumed that all the particles except the right-helicity neutrino are 
at a common temperature T and the chemical potentials of all the particles can be 
neglected. With Eq. (|2^) we are able to rewrite (0) in the more compact form 



{C2->2 + C'l^a) 



A. 



(23) 



where f™ = [exp{E+/T) + 

Let us introduce the total production rates per unit volume: 

T2^2{R) = J2 J dTiadmil+dUd{27T)^6^^\pa+Pb-p+-pd) 
scatt 

X ^iA^;^,+.ivr(i - c)(i - /r) • (24) 



dec 



Note that the expressions given above do not include the blocking factors (1 — /+). 
Also, the dependence of F/'s on R arise through the distribution functions of ambient 



particles. Comparing Eqs. (^) and (p^ one easily sees that 

ri{R) = /^^KIP+I'^)' (25) 

or equivalently 

From this expression we see that the quantities Cj, which are of the dimension of 
mass, are the production rates of z/+'s per unit volume per unit interval of (67r^) 
around |p_,_|. 

Eq. ( 126|) suggests also a general method for calculating C/'s. One first generates 
by a Monte Carlo (MC) method a sample of unweighted events, together with their 
common weight, such that the total integral over the initial and final state phase 
space yields Tj. Having this sample, the evaluation of the derivative in (pB]) essentially 
reduces to building the histogram the number of unweighted events vs. |p_|_|'^/(67r^). 
Although this method is very general (applicable both to scatterings and decays) and 
with a clear physical meaning, it yields a good accuracy only if one generates the initial 
sample of weighted events such that there are more weighted events in the regions of 
the phase space where the integrand in Tj is larger. This requires some knowledge of 
the behaviour of the integrand and an ability to build a MC event generator with the 
required distribution of events in the phase space. 

In the case of 2 — ^> 2 scatterings there exists also another method of calculation 
of C2_,2- This method follows directly from (^) and is based on the T-invariance of 
interactions allowing to rewrite 6*2^2 in a form where z/+ is an initial state particle. 
After this transformation the calculation of the corresponding integral proceeds in a 
usual way. In order to check our results we have determined 6*2^2 with both methods. 

3.3 Results 

The total rates C(|p+|,-R) = C2~>2{\p+\, R) + C'l-^sdp+l, -R) ^^e right-helicity tau 
neutrino production for a number of fixed momenta are shown in Fig. 2. One can 
clearly see the effect of the s-channel pole, indicating that the rate is dominated by 
fermionic s-channel processes. Small momentum states pass through the pole at high 
temperatures, large at low temperatures. Another feature is the increase of the total 
rate for very small momentum z/+'s at 0.2 GeV < T < 2 GeV. This enhancement is due 
to both scatterings {ud z^+t+, cs i^+t^) and the decays of the tau lepton. Since at 
these temperatures the thermally averaged center of mass energy of colliding particles 
is of the same order of magnitude as the mass of the tau lepton, the neutrinos (and r's) 
are preferentially produced in the small-momentum states, for which the probability 
of the creation of a right-helicity neutrino (spin-flip) is larger. For similar reasons the 



contribution of r-decays to the total rate can be as large as ~ 30%. For right-helicity 
muon neutrinos one expects the analogous mechanisms to be effective at temperatures 
below 100 MeV. 

Because different momentum states have different interaction rates, this causes a 
distortion of the momentum distribution of z/+'s relative to the equilibrium distribu- 
tion. This effect is demonstrated in Fig. 3. Here we have parametrised the distribution 
function of for fixed R, /+(|p+|), introducing a momentum- dependent effective tem- 
perature Teff(|p+|) through /+(|p+|) = [exp(E+/Teff(|p+|)) + 1]"^ In Fig. 3 the ratio 
Tefj/T for a tau neutrino of mass 10 keV is shown at T = 30, 3 and 0.3 GeV for two 
extreme initial conditions at a high temperature T = 100 GeV: (z) /+ = /_ = /™, 
i.e. complete equilibrium, and {ii) /+ = 0. One sees that only the lowest momentum 
right-helicity tau neutrinos of this mass will come into full equilibrium at the pole: 
Teff depends on the initial condition even at T = 0.3 GeV for |p+|/T > 1. The rise 
of Tgfj at the small-momentum end of the spectrum is because these states, interact- 
ing more strongly, are kept in good thermal contact with ambient matter even when 
higher momentum states are (already) decoupled. This kind of the distortion of the 
distribution function becomes more prominent as the mass of the neutrino increases. 
As mentioned before, the regime of free expansion reveals itself in the shape-preserving 
evolution of the distribution function. Another remarkable feature is the rise of T^s at 
higher momenta clearly seen for the curve corresponding to T = 30 GeV and the zero 
initial condition /_|_ = 0. This rise is there despite the fact the total production rate 
C(|p+|, -R) is a monotonically (at least for higher momenta) dereasing function of |p+|. 
The explanation can be found from Eq. (|23|) : if /+//™ is much smaller than unity, 
the destruction of z/+'s can be neglected. The regime of well-out-of-equilibrium pro- 
duction of Z/+ is effectively maintained only while both the interaction rate C(|p+|, R) 
and the ratio f+/ are sufficiently small. These conditions are more easily fulfilled 
by the higher momentum neutrinos at the beginning of their evolution from the zero 
initial condition. This out-of-equilibrium production results in the distribution func- 
tion decreasing with |p_|_| more slowly than the equilibrium distribution function, or, 
equivalently, an Tgfj rising with |p+|. In contrast with the above discussed rise of 
Teff at lower momenta, the later feature becomes more apparent as the mass of the 
neutrino decreses. We find that spectral distortion changes the energy desity of i/+'s 
(which were in equilibrium at temperatures above the weak interaction pole) at ~ 1 
MeV typically by only a few percent. The major difference between the actual energy 
density and the naive estimate pn ~ 0.044 is due to the relatively late decoupling 
caused by the pole. 

Once the initial condition is fixed, we may follow the evolution of the energy density 
of In Fig. 4 we show the evolution of the right-helicity neutrino energy density for 
m^^ = 1, 6 and 20 keV for two different initial conditions. If m^^ = 1 keV, spin-flip 



interactions are too weak to affect the evolution, and the resulting relic abundance at 
nucleosynthesis is, to a high accuracy, what one would naively expect. However, if 
rriu^ = 20 keV, spin-flip interactions at the pole are strong enough to equilibrate z/+ 
even if initially /+ = 0. 

The right-helicity neutrino energy density at about nucleosynthesis time in units 
of left-helicity neutrino energy density is shown in Fig. 5. We see that almost full 
equilibration of z/_|_ is obtained at some temperature below 100 GeV if m^^ > 10 keV. 
We have assumed here that there is no entropy production at QCD phase transition 
so that the only effect is the dilution of the right-helicity neutrino densities by the 
appropriate ratios of the effective degrees of freedom before and after the phase tran- 
sition. We have computed AN^ assuming also that below Tqcd all interactions can 
be ignored. While this is a very good approximation for smaller masses {m^^ < 30 
keV), it has been found p, ^ that for higher masses out-of-equilibrium scatterings 
and the decay 7r° J-'+T^+ (and the decay 7r+ fJ'~^^+ in the case of muon neutrinos) 
produce considerable amounts of z/+. Using the results of we estimate that the non- 
equilibrium reactions occuring below Tqcd give an additional contribution to AA^",^, to 
be added to our result shown in Fig. 5, which is approximately 0.003 for m^^ = 30 
keV, 0.012 for m^^ = 60 keV, 0.03 for m^^ = 100 keV and 0.15 for m^^ = 200 keV, 
adopting Tqcd = 150 MeV. Given the inaccuracy of the determination of N^, from the 
observational data, at present the pole effect cannot be used to derive an upper bound 
on the masses of Dirac neutrinos. 

We have assumed no other interactions than those present in the Standard Model. 
Any non-standard interaction above the electroweak scale would affect the initial en- 
ergy density of the right-handed neutrinos. However, according to our results, if the 
mass of neutrino is more than about 10 keV, nucleosynthesis is not sensitive to such 
interactions as their effect will be washed out by equilibration at the electroweak pole 
region. Because this equilibration is inefficient below neutrino masses in the 1 keV 
range, the cosmological mass limit of stable Dirac neutrinos is not modified. The 
increase is at most about 4% as expected naively, depending on the unknown initial 
density at high temperatures. 
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Table 1: Processes, up to crossings, which produce a right-hehcity tau neutrino 



fermionic 2^2 



bosonic 2 — > 2 



3 — body decays 



2 — body decays 



r i/j- ^ , {j = e, //) 



T 7 — > 



T — > V-\-dnUm , 

{dnUm — du, su, sc, dc) 
t dnT'^iy+ , 
{dn = b, s, d) 



b Umr'^u+ , 



(Um = C, U) 



Z — > l^+Vj- 
H — > V-^-Vr 



Figure captions 



Figure 1. Thermally averaged production rates per one right-helicity tau neutrino i^^ 
of mass 20 keV for the s-channel process ud —>■ z^+r"*", the t-channel process t~u v^d 
and the bosonic process r~7 — > u+W~. The Hubble expansion rate is also shown for 
comparision. 

Figure 2. The total production rates (dashed curves) C(|p+|,i?) of right-helicity 
tau neutrinos as functions of the temperature T. Prom top to bottom, the rates are 
given for p+/100 GeV = 0.33, 1, 2, 4, 6, 8 and 10 . For comparison the evolution of the 
Hubble parameter H is also given (solid curve). 

Figure 3. The ratio of the momentum-dependent effective temperature Tcs of the 
right-helicity tau neutrinos of mass 10 kcV to the plasma temperature T at T = 
30 GcV, 3 GcV and 0.3 GcV. In each case the solid curve corresponds to the equilibrium 
initial condition /+ = /™, and the dashed curve to the zero initial condition /+ = 
at T = 100 GeV. 

Figure 4. The evolution of p+ / of tau neutrinos as a function of the temperature 
T for rrii,^ — 20 keV (sohd curve), 6 keV (dot-dashed curve) and 1 keV (dashed curve). 
The evolution is given for two different initial conditions: /+ = /™ and /+ = at 
T = 100 GeV. 

Figure 5. The energy density of right-hehcity tau neutrinos u+ in units of the effective 
number of two-component neutrino species, AA^,^, as a function of the neutrino mass. 

ANi, is calculated at T = 3 MeV (solid curves), 1 MeV (dashed curves)and 0.3 MeV 
(dot-dashed curves). For rrii,^ < 20 keV the two sets of curves differ by initial conditions 
at T = 100 GeV. 
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